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This  paper  focuses  on  the  combination  of  wind  turbine  (WT),  photovoltaic  (PV),  fuel  cell  (FC)  and  ultra¬ 
capacitor  (UC)  systems  for  grid-independent  applications.  The  dynamic  behavior  of  the  proposed  hybrid 
system  is  tested  under  various  wind  speed,  solar  radiation  and  load  demand  conditions.  The  developed 
model  and  its  control  strategy  exhibit  excellent  performance  for  the  simulation  of  a  complete  day.  In  the 
simulation,  the  solar  radiation  and  power  demand  data  are  based  on  real  world  measurements,  while  the 
wind  speed  data  are  quasi-real  because  it  is  simulated  based  on  special  wind  speed  generation  algorithms. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Renewable  energy  sources  and  conversion  units  such  as  PVs, 
WTs,  and  FCs  have  been  found  to  be  promising  energy  sources 
toward  building  a  sustainable  and  environment  friendly  energy 
economy  in  the  next  decade.  However,  these  renewable  energy 
sources  suffer  from  some  deficiencies  when  used  as  stand-alone 
energy  sources.  The  power  generated  by  WT  and  PV  systems  is 
highly  dependent  on  weather  conditions.  Natural  variations  in 
wind  speed  and  sunlight  causes  power  fluctuations  in  WT  and  PV 
systems,  respectively.  In  addition,  it  is  difficult  to  store  the  power 
generated  by  a  PV  or  WT  system  for  future  use.  To  alleviate  these 
problems,  WT  and  PV  sources  can  be  integrated  with  other  alter¬ 
native  systems  using  hybrid  topologies.  A  hybrid  power  system 
consists  of  a  combination  of  two  or  more  energy  sources,  convert¬ 
ers  and/or  storage  devices.  Thus,  higher  efficiency  can  be  obtained 
by  making  the  best  use  of  their  features  while  overcoming  their 
limitations  [1-4]. 

Alternate  energy  conversion  systems  such  as  PV  panels  and  wind 
turbines  can  be  combined  with  FC  power  plants  to  satisfy  sustained 
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load  demands.  An  FC  power  plant  uses  hydrogen  and  oxygen  to  con¬ 
vert  chemical  energy  into  electrical  energy.  In  this  study,  among  the 
various  types  of  FC  systems,  a  Proton  Exchange  Membrane  (PEM)  FC 
power  plant  is  used  because  PEMFC  power  plants  have  been  found 
to  be  especially  suitable  for  hybrid  energy  systems  [5]  since  they 
have  low  operating  temperature  (80-1000  °C)  and  fast  start-up. 
Currently,  a  5  kW  commercial  PEMFC  power  plant  is  operational  in 
the  authors’  lab.  Moreover,  PEMFC  is  suitable  for  using  the  hydrogen 
produced  from  the  electrolysis  of  water  since  hydrogen  is  free  from 
any  carbon  monoxide  [1].  However,  assisting  an  FC  hybrid  power 
plant  with  a  parallel  UC  bank  makes  economic  sense  when  satisfy¬ 
ing  the  peak  power  demands  or  transient  conditions  since  the  UC 
bank  supplies  the  extra  power,  thereby  reducing  the  size  and  cost  of 
the  FC  system.  Ultra-capacitors  are  electrical  energy  storage  devices 
(a  few  Farads  to  several  thousand  Farads  per  cell)  with  high  power 
densities  when  compared  to  batteries  [6].  Ultra-capacitor  banks 
can  be  used  for  short-term  energy  storage  due  to  their  high  cycling 
efficiency,  convenience  for  charging/discharging,  and  additionally 
to  meet  the  instantaneous  load  ripples/spikes.  Besides,  overloading 
fuel  cell  systems  may  cause  gas  starvation  thus  decreasing  its  per¬ 
formance  and  lifetime.  As  another  feature,  the  load  tracking  delays 
and  mismatches  of  the  FCs  should  be  compensated  by  an  auxiliary 
system  such  as  a  battery  or  ultra-capacitor  bank. 

In  this  paper,  the  dynamic  model,  design  and  simulation  of  a  new 
WT/PV/FC/UC  autonomous  hybrid  system  is  proposed.  To  achieve 
autonomous  operation  and  energy  management,  appropriate 
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Fig.  1.  Overall  system  configuration. 


components  and  a  control  system  is  used  to  achieve  adequate  per¬ 
formance  for  the  overall  system.  Here,  the  control  system  has  been 
designed  to  draw  the  maximum  available  energy  from  the  WT  and 
PV  sources.  The  excess  power  generated  by  renewable  sources  is 
used  for  UC  charging  and  for  the  electrolyzer  to  produce  hydrogen. 
When  the  energy  generated  from  the  WT  and  PV  systems  becomes 
insufficient  with  respect  to  the  load  demand,  the  stored  hydrogen  is 
fed  to  the  FC  system.  The  tracking  mismatches  of  the  FC  system  and 
the  load  demand  exceeding  the  power  generated  by  WT/PV/FC,  is 
supplied  by  the  UC  bank.  Using  the  developed  model,  the  purpose  of 
testing  the  overall  performance  of  the  complete  system  is  achieved. 
Besides,  the  dynamic  models  of  each  component  are  integrated 
using  the  proposed  topology  with  good  performance.  The  model 
allows  to  design  and  size  the  hybrid  system  for  a  variety  of  load¬ 
ing  and  meteorological  conditions.  The  developed  model  provides 
the  integration  of  renewable  energy  systems  with  energy  man¬ 
agement  strategy  to  overcome  the  deficiencies  of  each  renewable 
system.  Modeling  and  simulations  are  performed  using  MATLAB®, 
Simulink®  and  SimPowerSystems®  [7-9]  software  packages  and 
test  results  of  the  autonomous  operation  and  energy  management 
of  the  system  show  that  this  technique  performs  successfully. 

2.  System  description  and  methodology 

The  overall  configuration  of  the  hybrid  power  generation  sys¬ 
tem  is  shown  in  Fig.  1.  The  PV/FC/UC  hybrid  power  system  model 
reported  in  Ref.  [10]  and  the  dynamic  PV  system  model  reported  in 
Ref.  [11  ]  are  enhanced  and  integrated  with  the  wind  turbine,  elec¬ 
trolyzer,  storage,  and  load  models.  The  primary  energy  sources  of 
the  system  consist  of  a  10  kW  WT  and  a  3  kW  (peak)  PV  array.  The 
excess  energy  with  respect  to  the  load  requirement  is  used  for  UC 
charging  and  hydrogen  production.  The  excess  power  demand  of 
the  user  is  supplied  by  the  FC  system  up  to  5  kW.  When  FC  system 
is  in  operation,  not  only  the  excess  demand  over  5  kW  but  also  the 
tracking  mismatches  and  delays  of  the  FC  system  are  compensated 
by  the  UC  bank. 

The  dc/dc  converters  depicted  in  Fig.  1  include  two  cascade  con¬ 
nected  converters.  The  bus  side  converters  are  used  for  bus  voltage 
regulation  and  the  source  side  converters  are  used  for  reference 
power  tracking. 


2.1  Wind  speed,  wind  turbine  and  generator  models 

The  Wind  Turbine  Blockset,  developed  in  Ref.  [12],  is  used  in 
this  study  to  realize  and  employ  the  wind  speed,  wind  turbine  and 
generator  models.  The  model  of  the  wind  speed  profile  is  based  on 
Kaimal  spectra  [12].  This  model  takes  into  account  the  rotational 
turbulences  and  the  tower  shadow.  The  built-in  band-limited  white 
noise  generator  is  used  from  Simulink  and  the  wind  speed  is  calcu¬ 
lated  as  an  averaging  of  the  fixed-point  wind  speed  over  the  whole 
rotor.  The  model  of  the  wind  turbine  is  based  on  the  torque  coef¬ 
ficient  Cq  using  a  look-up  table  as  a  function  of  A  and  9  for  the 
variable  pitch  wind  turbine.  The  torque  coefficient  Cq  is  used  to 
directly  determine  the  aerodynamic  turbine  torque,  given  by 

Tcot  =  0.5ttpRv2Cq.  (1) 

The  parameters  used  in  the  wind  turbine  model  are  as  follows: 

Cq  torque  coefficient 
R  blade  radius 

A  tip  speed  ratio 

p  air  density 

6  pitch  angle 

v  wind  speed 

An  active  stall  controller  is  also  employed  to  determine  the  opti¬ 
mal  pitch  angle  and  to  limit  the  electrical  output  power  to  the 
nominal  mechanical  power.  The  cut-in  and  cut-off  wind  speed  of 
the  turbine  is  set  to  be  3  m  s-1  and  15  m  s-1 ,  respectively.  The  out¬ 
put  torque  of  the  wind  turbine  is  used  to  drive  the  squirrel  cage 
induction  generator  shaft.  In  the  induction  generator  model,  the 
steady-state  voltage  equations  are  used.  Using  this  approach,  the 
induction  machine  can  be  completely  characterized  in  terms  of 
steady-state  values  of  the  stator  and  rotor  currents,  phase  angle, 
electromagnetic  torque  and  electrical  power.  The  electrical  power 
which  is  supplied  to  the  dc  bus  through  the  ac/dc  and  dc/dc  convert¬ 
ers  is  obtained  from  the  stator  windings  of  the  induction  generator. 
All  detailed  information  on  the  components  of  WT  system  such 
as  wind  model,  active  stall  controller,  aerodynamic  turbine  model, 
and  induction  generator  models  can  be  found  in  Ref.  [12].  The  WT 
system  configuration  is  given  in  Fig.  2. 
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Fig.  2.  WT  system  configuration. 


2.2.  PV  system  characteristics  and  model 

A  PV  system  consists  of  many  cells  connected  in  series  and  paral¬ 
lel  to  provide  the  desired  output  terminal  voltage  and  current,  and 
exhibits  a  nonlinear  I-V  characteristic  [11,13].  Masoum  et  al.  [11] 
introduced  a  model  for  Silicon  solar  PV  panel.  The  PV  cell  equivalent 
model,  which  represents  the  dynamic  nonlinear  I-V  characteristics 
of  the  PV  system  shown  in  Fig.  3,  is  used  within  the  hybrid  sys¬ 
tem  as  a  controlled  voltage  source  by  utilizing  it  in  the  Simulink 
and  SimPowerSystems  toolboxes  of  MATLAB.  The  subsystems  of  the 
complete  PV  model  are  described  in  details  in  Ref.  [10],  which  is  a 
previous  study  of  the  authors,  hence;  those  are  not  repeated  in  this 
study. 

The  parameters  used  in  the  mathematical  modeling  of  the  PV 
system  are  as  follows: 


a  ideality  or  completion  factor 

J0  PV  cell  reverse  saturation  current  [A] 

Jpv  PV  cell  output  current  [A] 

/sc  short-circuit  cell  current  [A] 

k  Boltzmann’s  constant  [J  K_1  ] 

Np  the  number  of  parallel  strings 

Ns  the  number  of  series  cells  per  string 
q  electron  charge  [C] 

Rs  series  resistance  of  PV  cell  [  £2  ] 

T  PV  cell  temperature  [K] 

VpV  terminal  voltage  for  PV  cell  [V] 


The  output  voltage  characteristic  of  the  PV  system  can  be 
expressed  as  in  Ref.  [11]  and  [13], 


VpV  = 


NsakT  , 

- In 

Q 


he  ~  /pv  +  Np 

Wo 


wfspv' 


(2) 


The  manufacturer’s  datasheet  provides  necessary  information 
for  most  of  the  parameters  ofEq.  (2)  [14].  The  current-voltage  char¬ 
acteristics  of  the  PV  array  (for  Np  =  30  and  Ns  =  112)  are  obtained 
according  to  the  different  short-circuit  current  values  using  Eq.  (2) 
as  shown  in  Fig.  4. 

Using  the  current-voltage  curves,  the  current-power  curves  can 
be  obtained  as  shown  in  Fig.  5  to  determine  the  maximum  power 
to  be  drawn  from  the  PV  array  under  various  short-circuit  current 
values. 

The  maximum  power  output  of  the  PV  array  varies  accord¬ 
ing  to  solar  radiation  or  load  current.  Therefore,  a  control  system 
is  needed  to  exploit  the  solar  array  more  effectively  as  an  elec¬ 
tric  power  source  by  building  a  maximum  power  point  tracker 
(MPPT)  [15,16].  Different  maximum  power  point  tracker  systems 
can  be  implemented.  In  this  case,  a  computational  approach  is  used 
involving  the  model  of  the  nonlinear  I-V  characteristics  of  the  solar 
panel.  Based  on  the  I-V  characteristics,  the  corresponding  maxi¬ 
mum  power  points  are  computed  for  different  load  conditions  as 
a  function  of  cell  short-circuit  currents  which  varies  as  a  function 
of  the  solar  radiation.  To  determine  the  relationship  between  the 
solar  radiation  and  short-circuit  current,  a  curve  fitting  based  com¬ 
putational  method  is  employed.  Then,  an  MPPT  system,  where  the 


sol  dat 


From 


Workspace 

Saturation 


1 93 


MosferGate  (  t 

l_out 

l_sc 

PV  Voltage 

- 

+ 

PV  System  Model 


solar  radiation 

Isc 

PV  Voltage 

Max.  Pow. 

w 

l_out 

MPPT  Gate 

V 

MPPT 


PV  Voltage 
l_out  Gate  Signal 
I  sc 


PV  Model 
Internal  Control 


dc/dc  Converter 
Controller 


Vdcjn 

Vdcjef 


dcBus 

Reference  Voltage 


0  3  3 

0  3  3 

Q.  Cl 

C  3 

MPPT 

Q.  Q. 

C  3 

O 

dc/dc 

o 

~o 

converter 

■a 

O 

O 

dc/dc 

Converter 


to 


dc  bus 


Fig.  3.  PV  system  with  MPPT,  controllers  and  converters. 
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l-V  Curves  of  the  Photovoltaic  Array  for  Different  I  Values 


Fig.  4.  Current-voltage  curves  for  different  short-circuit  levels. 


hardware  part  is  based  on  a  dc/dc  converter  is  built  to  draw  all  the 
maximum  available  power  from  the  PV  generator. 

The  voltage  based  PV  model  exhibits  the  l-V  characteristics  of 
the  PV  system  using  a  built-in  dc/dc  converter  as  shown  in  Fig.  3 
(PV  system  model).  The  pulse  width  of  the  switching  signal  is  deter¬ 
mined  using  reference  and  measured  voltages  to  represent  the 
characteristic  given  in  Eq.  (2).  Then,  the  pulse  width  command  and 
the  carrier  signal  are  compared  to  obtain  the  switching  PWM  sig¬ 
nals.  The  MPPT  system  calculates  the  maximum  available  power 
from  the  PV  system,  and  then  compares  the  calculated  power  with 
the  actual  drawn  power  for  generating  the  switching  signals  of  the 
dc/dc  converter  of  the  MPPT  system.  Another  dc/dc  converter  cas¬ 
caded  after  the  MPPT  converter  which  helps  to  keep  the  output 
voltage  equal  to  the  dc  bus  voltage. 

The  values  of  the  PV  system  parameters  are  listed  in  Table  1. 

2.3.  PEMFC  model 

The  PEMFC  model  described  in  Ref.  [17]  is  modified  using  the 
relationship  between  output  voltage  and  partial  pressure  of  hydro¬ 
gen,  oxygen  and  water.  The  FC  system  model  parameters  used  in 
this  model  are  as  follows: 


l-P  Curves  of  the  Photovoltaic  Array  for  Different  I  Values 


Fig.  5.  Current-power  curves  for  different  short-circuit  current  values. 


Table  1 

PV  system  model  parameters  [11] 


Specification 

Value 

The  number  of  series  cells  per  string  (Ns) 

30 

The  number  of  parallel  cells  per  strings  (Np) 

112 

Ideality  or  completion  factor  (a) 

1.9 

Boltzmann’s  constant  ( k ) 

1.3805e-23  [JR”1] 

PV  cell  temperature  (T) 

298  [K] 

Electron  charge  (q) 

1.6e-19C 

Short-circuit  cell  current  (representing  insulation  level)  (Jsc) 

2.926  [A] 

PV  cell  reverse  saturation  current  (Jo) 

0.00005  [A] 

Series  resistance  of  PV  cell  (Rs) 

0.0277  [Q] 

B,C  constants  to  simulate  the  activation  over  voltage  in  PEMFC 

system  and  [V] 

CV  conversion  factor  [kmol  of  hydrogen  per  kmol  of 
methane] 

E  Nernst  instantaneous  voltage  [V] 

E0  standard  no  load  voltage  [V] 

F  Faraday’s  constant  [C  (kmol)-1  ] 

JFC  FC  system  feedback  current  [A] 

/c2  PI  gain 

/Can  anode  valve  constant  [  y^kmolkg^tm  s)-1  ] 

I<h2  hydrogen  valve  molar  constant  [  y/kmoTkg(atm  s)-1  ] 
/Ch2o  water  valve  molar  constant  [  y/kmolkg(atm  s)_1  ] 

I<q2  oxygen  valve  molar  constant  [  y/kmol  kg(atm  s)_1  ] 

I<r  modeling  constant  [kmol  (s  A)-1  ] 

Mh2  molar  mass  of  hydrogen  (kg  kmol-1 ) 

N0  number  of  series  fuel  cells  in  the  stack 
Ph2  hydrogen  partial  pressure  [atm] 

Ph2o  water  partial  pressure  [atm] 

Po2  oxygen  partial  pressure  [atm] 

Qh2  molar  flow  of  hydrogen  [kmol  s-1  ] 
qo2  input  molar  flow  of  oxygen  [kmol  s-1  ] 

^methane  methane  flow  rate  [kmol  s-1  ] 
ql£2  hydrogen  input  flow  [kmol  s-1  ] 

q hydrogen  output  flow  [kmol  s-1  ] 
qrH2  hydrogen  flow  that  reacts  [kmol  s-1  ] 
qj^q  amount  of  hydrogen  flow  required  to  meet  the  load 
change  [kmols-1] 

R  universal  gas  constant  [(1  atm)  (kmol  K)-1  ] 

Rint  FC  internal  resistance  [£2] 

rH_o  the  hydrogen-oxygen  flow  ratio 

T  absolute  temperature  [K] 

U  utilization  rate 

Van  volume  of  the  anode  [m3  ] 

Vcen  dc  output  voltage  of  FC  system  [V] 

Th2  hydrogen  time  constant  [s] 

ro2  oxygen  time  constant  [s] 

Th2o  water  time  constant  [s] 

Pact  activation  over  voltage  [V] 

pohmic  ohmic  over  voltage  [V] 


The  relationship  between  the  molar  flow  of  any  gas  (hydrogen) 
through  the  valve  and  its  partial  pressure  inside  the  channel  can  be 
expressed  as  [18] 


qlh 

pH2 


(3) 


For  hydrogen  molar  flow,  there  are  three  significant  factors: 
hydrogen  input  flow,  hydrogen  output  flow,  and  hydrogen  flow  dur¬ 
ing  the  reaction  [18].  The  relationship  among  these  factors  can  be 
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expressed  as 


(4) 


According  to  the  basic  electrochemical  relationship  between  the 
hydrogen  flow  and  the  FC  system  current,  the  flow  rate  of  reacted 
hydrogen  is  given  by  [18] 


nr  _  NoNs/fC 
<?h2  -  2 F 


=  2  JCr/rc. 


(5) 


Using  Eqs.  (3)  and  (5),  and  applying  Laplace  transform,  the  hydrogen 
partial  pressure  can  be  obtained  in  the  s  domain  as  [19] 


P  h2 


yi<u2 

1  +  ?w2s 


(<? 


h2 


2/Cr/pc), 


(6) 


where 

Van 

H2  /Ch2JRT' 


(7) 


Similarly,  water  partial  pressure  and  oxygen  partial  pressure  can 
be  obtained.  The  polarization  curve  for  the  PEMFC  is  obtained  from 
the  sum  of  the  Nernst’s  voltage,  the  activation  over  voltage,  and 
the  ohmic  over  voltage.  Assuming  certain  temperature  and  oxygen 
concentration,  the  FC  output  voltage  may  be  expressed  as  [17,20] 


l/ceii  —  E  -\-  T)  act  +  4ohmio  (8) 

where 


Table  2 

FC  system  model  parameters 

Specification  Value 


0.04777  [A-1  ] 

0.0136  [V] 

2 

96484600  [Ckmol-1] 
3.37  [s] 


Activation  voltage  constant  ( B ) 
Activation  voltage  constant  (C) 
Conversion  factor  (CV) 

Faraday’s  constant  (F) 

Hydrogen  time  constant  (rH2 ) 
Hydrogen  valve  constant  (JCH2 ) 
Hydrogen-oxygen  flow  ratio  (rH-o) 
Kr  constant  =  No/4F 
Methane  reference  signal  (Qmethref) 
No  load  voltage  (Eo) 

Number  of  cells  (No) 

Number  of  stacks  (Ns) 

Oxygen  time  constant  (ro2 ) 

Oxygen  valve  constant  ( ko2 ) 

FC  system  internal  resistance  ( Rint ) 
FC  absolute  temperature  ( T ) 
Universal  gas  constant  ( R ) 
Utilization  factor  (U) 

Water  time  constant  (rH2o) 

Water  valve  constant  (1<h2o) 


4.22  x  10-5  [kmol  (s  atm)-1  ] 
1.168 

2.2802  xlO-7  [kmol  (s  A)-1] 
0.000015  [kmols-1) 

0.8  [V] 

88 

1 

6.74  [s] 

2.11  x  10-5  [kmol  (s  atm)-1] 
No  x  0.00303  [Q] 

343  [K] 

8314.47  [J(kmolK)-1] 

0.8 

18.418  [s] 

7.716  x  10-6  [kmol  (s  atm)-1  ] 


Eye-total  the  total  UC  system  resistance  [£2] 

T  time  [s] 

V{t)  the  voltage  at  time  t  [V] 

Vi  the  initial  voltage  before  discharging  starts  [V] 

Vf  the  final  voltage  after  discharging  ends  [V] 


Pact  =  -B  In (C/pc) 

0) 

and 

*7ohmic  —  E  /pc- 

(10) 

may  be  expressed  as  [17] 

(11) 

according  to  the  power 
the  high-pressure  hydro¬ 
gen  tanks.  During  operational  conditions,  to  control  hydrogen  flow 
rate  according  to  the  output  power  of  the  FC  system,  a  feedback 
control  strategy  is  utilized.  To  achieve  this  feedback  control,  FC  cur¬ 
rent  from  the  output  is  taken  back  to  the  input  while  converting  the 
hydrogen  into  molar  form  [17,21  ].  The  amount  of  hydrogen  available 
from  the  hydrogen  tank  is  given  by 


Now,  the  Nernst’s  instantaneous  voltage 


E  =  N0 


RT  , 

Eo+2 Fl0g 


P  h2o 


The  FC  system  consumes  hydrogen 
demand.  The  hydrogen  is  obtained  from 


„req  _  N0NS/FC 
Qh2  ~  2FU 


(12) 


Depending  on  the  FC  system  configuration,  and  the  flow  of 
hydrogen  and  oxygen,  the  FC  system  produces  the  dc  output 
voltage.  The  hydrogen-oxygen  flow  ratio  rH_o  in  the  FC  system 
determines  the  oxygen  flow  rate  [17]. 

The  values  of  the  FC  system  parameters  are  listed  in  Table  2. 


2.4.  Ultra-capacitor  model 


Fig.  6  shows  the  classical  equivalent  circuit  of  the  UC  unit.  The 
model  consists  of  a  capacitance  (C),  an  equivalent  series  resis¬ 
tance  (ESR,  R )  representing  the  charging  and  discharging  resistance, 
and  an  equivalent  parallel  resistance  (EPR)  representing  the  self¬ 
discharging  losses  [21,22].  The  EPR  models  the  leakage  effects  and 
only  impacts  long-term  energy  storage  performance  of  the  UC 
[23,24].  The  voltage  state  of  an  UC  (i.e.,  RC  circuit)  with  a  capac¬ 
itance  C  draining  charge  into  a  resistance  R  may  be  expressed  as 
[25] 

V(t)  =  V iexp^-jU).  (13) 

The  RC  time  constant  determines  the  effective  period  of  the  charg¬ 
ing  and  discharging  processes  for  some  initial  voltage  on  the 
capacitor  [25].  The  state  of  charge  (SOC)  of  an  UC  bank  is  the  avail¬ 
able  energy  capacity  expressed  as  a  percentage  of  its  rated  energy 
capacity  which  is  related  to  its  voltage.  The  amount  of  energy  drawn 
from  the  UC  bank  is  directly  proportional  to  the  capacitance  and  the 
change  in  terminal  voltage  [23,6],  which  is  expressed  as 

£uc=lc(V2-Vf2).  (14) 

The  effective  specific  energy  for  a  prescribed  load  can  be  sup¬ 
plied  by  various  UC  bank  configurations.  In  practical  applications, 
the  required  amount  of  terminal  voltage  and  energy  or  the  capaci¬ 
tance  of  UC  storage  system  can  be  built  using  multiple  UCs  in  series 


The  parameters  used  in  the  mathematical  modeling  of  the  UC 
are  as  follows: 

C  capacitance  [F] 

Cue-total  the  total  UC  system  capacitance  [F] 

EPR  equivalent  parallel  resistance  [£2] 

ESR,R  equivalent  series  internal  resistance  [£2] 

E\jc  the  amount  of  energy  drawn  from  the  UC  bank  [Ws] 
ns  the  number  of  capacitors  connected  in  series 
np  the  number  of  series  strings  in  parallel 


Fig.  6.  Classical  equivalent  model  for  the  UC  unit. 
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Table  3 

Maxwell  Boostcap  BMOD2600-48EA  ulta-capacitor  specifications 


Specification 

Value 

Capacitance 

144  [F] 

Voltage 

48.6  [V] 

ESR,  dc 

11  [m£2] 

ESR  @ 1  kHz 

6.8  [m£2] 

Thermal  resistance 

0.27  pCW-1] 

Short-circuit  current 

4800  [A] 

Energy  density  (max) 

3.40  [Whkg-1  ] 

Power  density 

2600  [Wkg-1] 

Volume 

13.4  [L] 

Mass 

13.5  [kg] 

and  parallel.  The  terminal  voltage  determines  the  number  of  capaci¬ 
tors  which  must  be  connected  in  series  to  form  a  bank  and  the  total 
capacitance  determines  the  number  of  capacitors  which  must  be 
connected  in  parallel  in  the  bank.  The  total  system  resistance  and 
the  total  system  capacitance  of  the  UC  bank  can  be  calculated  as 
[21,23] 


2.5.  Electrolyzer  and  storage  models 

Water  can  be  decomposed  into  its  elementary  components  by 
passing  electric  current  between  two  electrodes  separated  by  an 
aqueous  electrolyte  [27-29].  The  electrochemical  reaction  of  water 
electrolysis  is  given  by 

H20(iiquid)  +  electrical  energy  =>  H2  +  ^02.  (17) 

The  Electrolyzer  model  involves  the  following  parameters: 

F  Faraday  constant  [  C  kmol- 1  ] 

ze  electrolyzer  current  [A] 

nc  the  number  of  electrolyzer  cells  in  series 

i) F  Faraday  efficiency 

Hh2  produced  hydrogen  moles  per  second  [mol  s-1  ] 

According  to  Faraday’s  law,  hydrogen  production  rate  of  an  elec¬ 
trolyzer  cell  is  directly  proportional  to  the  electrical  current  in  the 
equivalent  electrolyzer  circuit  [27],  which  can  be  expressed  as 


^UC-total  =  ns 


ESR 


Cue-total  = 


(15) 

(16) 


In  this  research,  Maxwell  Technologies  Boostcap®  BMOD2600- 
48 EA  type  ultra-capacitor  bank  is  used  which  is  currently 
operational  in  the  authors’  lab.  The  reasons  for  selecting  this 
module  are  the  ultra  low  internal  resistance,  low  RC  time  con¬ 
stant  and  high  power  performance.  In  this  study,  4  parallel  UC 
modules  are  used,  where  each  module  includes  18  UC  units  con¬ 
nected  in  series.  The  specifications  of  the  UC  bank  [26]  is  given  in 
Table  3. 


Hh2 


dr  nde 
2  F 


(18) 


The  ratio  between  the  actual  and  the  theoretical  maximum  amount 
of  hydrogen  produced  in  the  electrolyzer  is  known  as  Faraday  effi¬ 
ciency.  Assuming  that  the  working  temperature  of  the  electrolyzer 
is  40  °C,  Faraday  efficiency  is  expressed  by  Refs.  [27,29] 

rj  F  =  96.5  e(0-09/ie_75-5/'i).  (19) 


According  to  Eqs.  (18)  and  ( 19),  a  simple  electrolyzer  model  is  devel¬ 
oped  using  Simulink. 

The  amount  of  hydrogen  required  by  the  PEMFC  is  sent  directly 
from  the  electrolyzer  system  according  to  the  relationship  between 
the  output  power  and  the  hydrogen  requirement  of  the  PEMFC 
system.  The  remaining  amount  of  hydrogen,  i.e.,  the  difference 


Power  generated  by  Power  generated  by  Power  demand 


Fig.  7.  The  general  decision  algorithm  used  in  the  main  controller. 
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Fig.  10.  Power  produced  by  the  wind  generator. 


Fig.  8.  Power  demand  of  the  micro-grid  neighborhood  consisting  of  3  residential 
units. 


between  the  produced  and  consumed  hydrogen  is  sent  to  the  stor¬ 
age  tank. 

One  of  the  hydrogen  storage  techniques  is  physical  hydrogen 
storage,  which  involves  using  tanks  to  store  either  compressed 
hydrogen  gas  or  liquid  hydrogen.  The  hydrogen  storage  model  is 
based  on  Eq.  (20)  and  it  directly  calculates  the  tank  pressure  using 
the  ratio  of  hydrogen  flow  to  the  tank.  The  produced  hydrogen  is 
stored  in  the  tank,  whose  system  dynamics  can  be  expressed  as  [30] 


^b  ^bi 


^Nu2RTb 


(20) 


The  parameters  used  in  the  hydrogen  storage  system  are  listed 
below: 


Mh2  Molar  mass  of  hydrogen  [kg  kmol-1  ] 

Nh2  hydrogen  moles  per  second  delivered  to  the  storage  tank 
[kmols-1  ] 

Pb  pressure  of  tank  [Pa] 

Pbi  Initial  pressure  of  the  storage  tank  [Pa] 

R  universal  (Rydberg)  gas  constant  [J  (kmol  I<)-1  ] 

Tb  operating  temperature  [K] 

Vb  volume  of  the  tank  [m3  ] 

z  compressibility  factor  as  a  function  of  pressure 


The  compression  dynamics  and  the  compression  energy 
requirements  are  not  considered  in  our  calculations.  All  auxil¬ 
iary  power  requirements  such  as  pumps,  valves  and  compression 
motors  are  ignored  in  the  dynamic  model  since  their  investigation 
is  out  of  the  paper’s  scope. 

2.6.  Energy  management  and  control  strategy 

A  main  controller  is  employed  for  the  entire  energy  measure¬ 
ment  and  management  purposes.  The  main  controller  provides 
the  autonomous  operation  with  the  required  measurements,  deci¬ 
sions  and  controls  by  collecting  data  through  the  sensors  and 
producing  the  commands  for  the  power  converters  connected  to 
the  inputs/outputs  of  the  components  used  in  the  hybrid  system. 
The  operating  strategies  used  in  the  main  controller  are  as  fol¬ 
lows: 

-  The  use  of  power  generated  by  the  PV  and  WT  systems  has  priority 
in  satisfying  power  demand  over  that  provided  by  the  FC  system 
or  by  the  UC  bank. 

-  If  the  total  power  generated  by  the  PV  and  WT  systems  is  higher 
than  the  demand,  the  additional  power  will  be  used  to  charge  the 
UC  bank. 

-  After  charging  the  UC  bank,  the  remaining  power  is  used  for  the 
electrolyzer  through  a  power  converter  to  generate  hydrogen. 


Fig.  9.  Wind  speed  profile. 


Fig.  11.  Global  hourly  solar  radiation. 
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lsc  -  P  Relationship  of  the  PV  Cell 


Fig.  12.  Short-circuit  current  versus  power  variation. 

-  The  electrolyzer  is  kept  in  operation  as  long  as  available  power 
exists  and  the  tank  pressure  does  not  exceed  the  prescribed  value 
of  15,000  Pa.  The  electrolyzer  maximum  power  is  also  set  to 
10  kW.  The  additional  power  generated  by  the  WT  and  PV  systems 
can  be  transferred  to  the  electrolyzer  up  to  the  10  kW  limit. 

-  Dumping  power  is  necessary  if  the  UC  bank  is  fully  charged  and 
the  electrolyzer  power  is  at  10  kW  limit. 

-  If  the  total  electric  power  generated  by  the  PV  and  WT  systems 
is  less  than  the  demand,  power  will  be  supplied  from  the  FC  sys¬ 
tem  up  to  the  maximum  5  kW  limit.  If  the  load  demand  exceeds 
the  power  generated  by  WT/PV/FC  combination,  the  difference  is 
supplied  by  the  UC  bank. 

-  The  UC  bank  not  only  provides  the  excess  power  demand  over 
5  kW,  but  also  compensates  the  tracking  mismatches  and  delays 
of  the  FC  system  which  has  relatively  slow  response  time. 

-  The  UC  bank  maximum  voltage  is  set  to  be  its  nominal  volt¬ 
age  of  48.6  V  and  its  minimum  voltage  is  set  to  be  25%  of 
full  state  of  charge  (SOC)  to  avoid  overcharging  and  deep 
discharging. 

To  implement  the  aforementioned  power  management  strategy, 
appropriate  power  control  systems  are  used  at  relevant  points  and 
the  total  system  components  are  integrated.  The  block  diagram  of 
the  main  controller  is  illustrated  in  Fig.  7. 


Fig.  14.  The  difference  between  the  load  power  and  generated  power  (PV  + WT). 

3.  Results  and  discussions 

In  the  simulation  process,  the  aim  is  to  observe  the  proposed 
system’s  behavior  over  a  long  period  of  time  including  day  and  night 
cases  as  well  as  low  or  no  wind  conditions.  The  load  profile,  wind 
speed  variation  and  solar  radiation  profiles  are  all  used  to  test  the 
performance  of  the  proposed  hybrid  system  for  a  typical  day  of  the 
neighborhood. 

The  5  kW  PEMFC  system  currently  operated  in  our  laboratory 
feeds  a  500  ft2  house.  Energy  consumption  varies  from  house  to 
house,  but  to  simulate  a  real-world  scenario,  load  profiles  are  also 
obtained  from  a  neighborhood  comprising  of  three  houses,  where 
each  house  is  approximately  2500  ft2  and  includes  all  household 
electrical  devices.  The  load  profile  of  the  micro-grid  neighborhood 
is  illustrated  in  Fig.  8. 

The  wind  speed  profile  is  generated  by  the  wind  speed  model 
described  in  Section  2.1  and  the  hourly  wind  speed  profile  is 
depicted  in  Fig.  9,  which  yields  the  power  generated  by  WT  as 
shown  in  Fig.  10. 

The  hourly  averaged  global  solar  radiation  on  a  horizontal  sur¬ 
face  (W  m-2)  is  taken  from  Ref.  [31  ],  where  the  record  coordinates 
are:  Latitude  41  °10'  North  and  Longitude  40° 20'  East.  The  exper¬ 
imental  solar  radiation  data  which  is  recorded  by  the  regional 
meteorological  station  in  the  month  of  June  is  given  in  Fig.  11. 


Fig.  13.  Hourly  power  drawn  from  the  solar  system. 


Fig.  15.  Power  produced  by  the  FC  system. 
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Time  [hour] 


Fig.  16.  Tracking  mismatches  of  the  FC  system. 


Fig.  18.  Amount  of  hydrogen  produced  by  the  electrolyzer. 


Using  the  nonlinear  curve  fitting  based  computation,  the  rela¬ 
tionship  between  the  solar  radiation  and  the  short-circuit  current 
is  obtained.  Then  using  the  maximum  power  points  for  different 
short-circuit  values  given  in  Fig.  5,  the  short-circuit  current  and  the 
maximum  power  relationship  are  obtained  as  shown  in  Fig.  12.  It 
is  evident  from  Fig.  12  that  the  short-circuit  current  varies  almost 
linearly  with  the  solar  radiation. 

The  maximum  power  values  are  used  in  the  MPPT  dc/dc  con¬ 
verter  as  power  commands  and  the  maximum  power  drawn  from 
the  PV  system  to  the  dc  bus  is  obtained  as  illustrated  in  Fig.  13.  From 
Fig.  13,  we  observe  that  in  the  early  morning  before  4:08  and  in 
the  evening  after  19:38,  the  solar  power  is  unavailable  due  to  non¬ 
existence  of  solar  radiation.  The  maximum  power  from  PV  system 
can  be  transferred  to  the  dc  bus  at  about  12.02  in  the  noon. 

As  mentioned  earlier,  the  use  of  power  generated  by  the  PV  and 
WT  systems  has  priority  in  satisfying  load  demand.  Thus,  the  total 
power  generated  by  the  WT  and  PV  systems  is  transferred  to  the  dc 
bus  directly.  The  difference  between  the  power  generated  by  WT 
and  PV  systems  and  the  user  power  demand  is  demonstrated  in 
Fig.  14.  The  positive  regions  of  Fig.  14  are  the  additional  power  gen¬ 
erated  by  PV  and  WT  sources  which  is  used  for  the  UC  charging  and 
hydrogen  production  within  the  electrolyzer  system.  Furthermore, 
the  negative  regions  are  the  exceeded  power  demands  which  are 
met  by  the  FC  system  and  UC  bank  discharge. 


The  FC  system  power  command  is  obtained  according  to  the 
negative  regions  of  Fig.  14  (exceeded  user  demands)  up  to  the  5  kW 
limit.  The  power  produced  by  FC  system  is  given  in  Fig.  15,  which 
slightly  differs  from  the  power  command  due  to  the  time  constants 
and  relatively  slow  response  characteristic  of  the  FC  system.  FIow- 
ever,  since  the  required  hydrogen  is  not  processed  by  a  reformer 
and  provided  directly  from  the  hydrogen  tank,  the  reformation  time 
delays  are  eliminated. 

The  tracking  mismatches  of  the  FC  system  are  given  in  Fig.  16 
which  is  supposed  to  be  supplied  by  the  UC  bank.  The  power  pro¬ 
duced  by  the  FC  system  shown  in  Fig.  15  corresponds  to  the  amount 
of  hydrogen  flow  to  the  FC  system  as  depicted  in  Fig.  17.  The  amount 
of  power  transferred  to  the  electrolyzer  corresponds  to  the  hydro¬ 
gen  flow  from  the  electrolyzer  to  the  storage  tank  as  illustrated  in 
Fig.  18. 

The  hydrogen  produced  by  the  electrolyzer  and  consumed  by  the 
FC  system  causes  the  pressure  of  the  storage  tank  to  vary  as  shown 
in  Fig.  19.  The  pressure  of  the  storage  tank  increases  when  the  elec¬ 
trolyzer  produces  hydrogen,  and  the  tank  pressure  decreases  when 
the  FC  system  uses  the  stored  hydrogen. 

The  power  supplied  by  the  UC  bank  is  shown  in  Fig.  20  which 
includes  the  exceeded  load  demand  and  the  tracking  mismatches 
(Fig.  16)  of  the  FC  system.  Due  to  the  scaling  differences  of 
Figs.  16  and  20,  the  low  mismatches  cannot  be  seen  in  Fig.  20.  These 


xio"4 


Fig.  17.  Amount  of  hydrogen  required  by  the  FC  system.  Fig.  19.  The  pressure  variation  of  the  hydrogen  storage  tank. 
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Fig.  20.  Power  supplied  by  UC  bank  including  tracking  delays  and  all  mismatches. 


Fig.  21.  UC  bank  charge  and  discharge  power. 


mismatches  can  be  observed  by  zooming  into  the  areas  of  interest 
of  Fig.  20. 

The  charge  and  discharge  power  of  the  UC  bank  is  recorded  as 
shown  in  Fig.  21 ,  where  positive  power  region  represents  the  power 
released  by  the  UC  bank  and  negative  power  region  represents  the 
power  captured  by  the  UC  bank  through  the  bi-directional  dc/dc 
converter. 


Fig.  22.  UC  bank  terminal  voltage. 


The  charging  and  discharging  power  of  the  UC  bank  corresponds 
to  the  terminal  voltage  variation  as  demonstrated  in  Fig.  22.  From 
Fig.  22,  it  is  obvious  that  the  UC  bank  terminal  voltage  varies  with 
the  charge  and  discharge  state  of  the  UC  bank  within  the  limited 
values. 

4.  Conclusions 

A  novel  WT/PV/FC/UC  hybrid  power  system  is  designed  and 
modeled  for  grid-independent  applications.  A  detailed  simulation 
model  has  been  developed  which  allows  designing  and  analyzing 
any  WT/PV/FC/UC  hybrid  system  with  various  power  levels  and 
parameters.  The  dynamic  behavior  of  the  hybrid  system  is  tested 
under  various  wind  speed,  solar  radiation  and  load  demand  con¬ 
ditions.  The  solar  radiation  and  power  demand  data  are  based  on 
real  world  records  and  the  wind  speed  is  artificially  generated  in 
the  computer  environment  using  a  suitable  algorithm. 

The  main  contribution  of  this  work  is  the  hybridization  of  the 
renewable  energy  sources  with  FC  systems  using  long  and  short¬ 
term  storage  strategies  with  appropriate  power  controllers  and 
control  strategy  to  build  an  autonomous  system.  The  developed 
system  and  its  control  strategy  exhibit  excellent  performance  for 
the  simulation  of  a  complete  day  or  longer  periods  of  time.  The 
proposed  system  can  be  used  for  non-interconnected  remote  areas 
or  isolated  cogeneration  power  applications  as  a  reliable  and  tech¬ 
nically  feasible  system. 
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